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Using an electron-positron collision data sample of 2.93 fb `! collected at a center-of-mass energy 
of J/s = 3.773 GeV with the BESIII detector, we present the first search for the radiative leptonic 
decay D* — ye*ve. The analysis is performed with a double tag method. We do not observe a 
significant Dt — e^ ue signal, and obtain an upper limit on the branching fraction of Dt > ye v. 
decay with the energy of radiative photon larger than 10 MeV of 3.0 x 10 ` at the 90% confidence 


level. 


PACS numbers: 13.20.Fc, 12.39.St 


I. INTRODUCTION 


In contrast to the purely leptonic decay, the radiative 
leptonic decay of the charged charmed meson, Dt — 
ye ue, is not subject to the helicity suppression rule due 
to the presence of a radiative photon. With no final- 
state hadron, treatment of the non-perturbative strong 
interaction effects in theoretical calculations is relatively 
simple. 

The radiative leptonic decays of heavy mesons have 
been studied with various models [144. Within the per- 
turbative quantum chromodynamics (pQCD) approach, 
the branching fraction of Dt + yeT ue decay is predicted 
to be of order 1074 {i}. Much smaller branching fractions, 
of order 1079, are obtained within the light front quark 
model B and the non-relativistic constituent quark mod- 
el E In Ref. A, the long-distance contribution is con- 
sidered via the vector meson dominance (VMD) model 
and it is found that the decay rate may be enhanced 
significantly. To deal with non-perturbative effects, it is 
important to separate the hard and soft physics, typi- 
cally with an approach known as factorization. Many 
approaches to factorization of the radiative leptonic de- 
cays of heavy mesons have been proposed [5410]. In re- 
cent papers miei factorization is extended to consider 
the first-order corrections in the strong coupling constant 


as and the heavy quark mass; the branching fraction of 
D+ + *ye*v, decay is predicted to be of order 1075. 

In this paper, we present the first search for the deca; 
D+ + *e*v,, based on a data sample of 2.93 fb ^! TEIE 
collected with the BESIII detector at a center-of-mass 
energy ys = 3.773 GeV. No obvious signal is observed, 
and an upper limit on the branching fraction of Dt > 
ye v, decay is set at the 90% confidence level (C.L.). In 
this paper, charge conjugate modes are always implied. 


II. THE BESIII DETECTOR AND DATA SET 


'The BESIII detector is a general purpose spectrometer 
with a geometrical acceptance of 93% of 47. It consists 
of a main drift chamber (MDC) for measuring the mo- 
mentum and specific ionization of charged particles in 
a 1T solenoidal magnetic field, a time of flight (TOF) 
system to perform particle identification, and a CsI(TI) 
electromagnetic calorimeter (EMC) for measurement of 
deposited shower energies. These components are sur- 
rounded by a multi-layer resistive plate counter system, 
which is designed to identify the muons. A detailed de- 
scription of the BESIII detector can be found in Ref. [13]. 

High-statistics Monte Carlo (MC) simulated data sam- 
ples are used to determine the detection efficiency and 
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to estimate potential background contamination. A 
GEANT4-based MC simulation program is used to 
simulate the interactions of particles in the spectrom- 
eter and the detector response. For the production of 
v (3770), KKMC [17] is used; it includes the effects of beam 
energy spread and initial state radiation (ISR). The 
known decay modes are generated using EVTGEN [19] 
according to branching fractions from the Particle Data 
Group (PDG) [20], and the remaining unknown decay 
modes are simulated by LUNDCHARM [1]. Final-state 
radiation (FSR) of charged tracks is incorporated with 
PHOTOS [22]. In modeling the signal events, the approach 
of Ref. is adopted, where first-order effects in the 
strong coupling constant a, and the heavy quark mass 
are considered. The minimum energy of the radiative 
photon is set at 10 MeV to aoud Pus infrared divergence 
for soft photons. For Dt — n?etve decay, which is an 
important background, an exclusive MC sample is gen- 
erated by adopting the associated form-factor model and 
parameters in Ref. [23]. 


HI. Dt > yetu. DATA ANALYSIS 


The analysis uses a double-tag (DT) technique 
which exploits the exclusive DD final states produced 
near threshold in e*e^ experiments. This technique 
allows one to measure absolute decay branching frac- 
tions of D* mesons independent of any direct knowl- 
edge of the total number of D*D~ events. In this 
analysis, the D^ candidates, so-called single-tag (ST) 
events, are reconstructed through six specific Pon 
decay modes K*m - , KTm m m^, Kom, E T? 
K}łn*r r and K*K mw . The signal Dt + yetv, is 
then searched for among the remaining tracks and show- 
ers recoiling against the ST D^ candidates; such signal 
candidate events are denoted as double-tag (DT) events. 
The absolute branching fraction, B(D* — yetve), can 
be obtained from the ratio of the DT yields and the ST 
yields, 


Npr 
> Nérgbr/er 


where Npr is the sum of signals yields for all tag modes 
and Ni, £54 and ei, are the ST yields and the detec- 
tion efficiencies of DT and ST for ST mode i, respectively. 
With this approach, the systematic uncertainties in the 
ST selection reconstruction are largely canceled in the 
branching fraction measurement. 


B(D* + yetve) = (1) 


A. Single-Tag event selection and yields 


For each charged track, we require the polar angle 0 
in the MDC to satisfy | cos0| < 0.93 and the point of 
the closest approach to the interaction point (IP) of the 
e*e~ beams to be within 1 cm in the plane perpendicular 


to the beam (V,) and within +10 cm along the beam axis 
(Vz). Particle identification (PID) for charged tracks is 
accomplished by combining the information on the mea- 
sured ionization energy loss (dE/dx) in the MDC and 
the flight time in the TOF into a PID likelihood, L(h), 
for each hadron hypothesis h = K or h = m. The 
7 (K) candidates are required to satisfy L(t) > L(K) 
(L(K) > £(n)). 

The Kg candidates are reconstructed from combina- 
tions of two tracks with opposite charge which satis- 
fy |cos0| < 0.93 and |V;| < 20cm, but with no V, 
and no PID requirements. The KÌ candidates must 
have an invariant mass in the range 0.487 < Mr+r- < 
0.511 GeV/c?, corresponding to three times our mass res- 
olution. To reject combinatorial background, we further 
require the decay length of K2 candidates, the distance 
between the IP and the reconstructed secondary decay 
vertex provided by a vertex fit algorithm, to be larger 
than two standard deviations. The momenta of ntr 
pairs after the vertex fit are used in subsequent analysis. 

'Those showers deposited in the EMC not associated 
with charged tracks are identified as photon candidates. 
'The energy deposited in the nearby TOF counters is in- 
cluded to improve energy resolution and detection effi- 
ciency. The minimum deposited energy is required to be 
greater than 25 MeV in the barrel region (| cos 0| < 0.80), 
or 50 MeV in the end caps regions (0.84 < | cos 0| < 0.92). 
'The shower time is required to be within 700 ns after the 
event start time to suppress electronic noise and show- 
ers unrelated to the collision event. The z? candidates 
are reconstructed from pairs of photons with invariant 
mass satisfying 0.115 < My, < 0.150 GeV /c?; those with 
both photons in the EMC end caps are rejected because 
of poorer resolution. The photon pairs of z? candidates 
are subject to a one-constraint (1C) kinematic fit which 
constrains their mass to the nominal z? mass [2d]; the 
updated momenta are used in subsequent analysis. 

The ST D signals are discriminated from back- 
grounds based on two kinematic variables, the energy dif- 
ference, AE, and the beam-constrained mass, Mpc (en- 
compassing energy and momentum conservation) which 
are defined as: 


AE = Est = beam; (2) 


Mpc = \/ Best m [ps |? Je (3) 


where psr and Esr are the total momentum and energy 
of the ST D^ candidate in the rest frame of the ete 
system, respectively, and Epeam is the beam energy. The 
ST signals peak around zero in the AF distribution and 
around the nominal D^ mass in the Mpc distribu- 
tion. 

For each ST mode, the D~ candidates are reconstruct- 
ed from all possible combinations of final-state particles, 
and are required to have AF within the regions listed in 
Table lll these are final-state dependent and determined 
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from data. If multiple candidates are found, only the 
one with the smallest |A E| is selected. To extract the 
ST signal yields, we perform extended unbinned maxi- 
mum likelihood fits to the Mgc distributions, as shown 
in Fig. [1| In the fits, signal shapes derived from the sig- 
nal MC events are convoluted with a Gaussian function; 
the free mean and width of this Gaussian compensate for 
imperfections in the beam energy calibration and differ- 
ences in the detector resolution between data and MC 
simulation, respectively. The combinatorial background 
is modeled by a smooth ARGUS function (25). The signal 
yields and the corresponding detection efficiencies in the 
region 1.8628 < Mpc < 1.8788 GeV/c? are summarized 
in Table [| A study of the inclusive DD MC samples, 
in which both D mesons decay inclusively, indicates that 
there are no significant backgrounds which peak in Mpc. 


TABLE I. Summary of the AE requirements, ST yields Når 
in data and detection efficiencies &&,. The efficiencies do not 
include the branching fractions of K} — «^«^ and 7° > yy. 
All uncertainties are statistical only. 


Tag Mode AE (MeV) Ni ei (96) 
Kin a [-27,25|] 801498 X 940 51.57 + 0.02 
Knn r [—62,34] 242092 + 699 24.37 + 0.02 
Kon [-25,25] 98132 + 328 54.03 + 0.06 
Kea a? [-73,41] 213976 + 641 26.17 + 0.02 
Knn m [—33,30] | 127463 + 415 32.46 + 0.04 
K+ K-14 [-23,20] 70701 + 343 


41.83 + 0.06 
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FIG. 1. (color online) The Msc distributions for the six tag 
modes. Dots with error bars are data. The blue solid lines 
show the overall fit curves and the red dashed lines are for 
the background contributions. 


B. Double-Tag event selection and yields 


We search for the signal Dt — yetve in the remain- 
ing charged tracks and showers recoiling against the ST 
D~ candidates. Exactly one good remaining charged 
track is required, with charge opposite to that of the 
ST D-. The track must be identified as an electron by 
combining the information from dE/dz, TOF, and the 
EMC. The PID £ is required to satisfy C(e) > 0 and 
L(e)/(L(e) + L(m)+L(K)) > 0.8. There must be at least 
one remaining photon to be selected as the candidate ra- 
diative photon; in the case of multiple candidates, the 
highest energy photon is used. However, we reject events 
in which any pair of photons satisfies x? < 20 in the 7° 1C 
kinematic fit. To improve the degraded momentum res- 
olution of the electron due to FSR and bremsstrahlung, 
the energy of neighboring photons, presumably due to 
FSR, is added back to electron candidates. Specifically, 
photons with energy greater than 50 MeV and within a 
cone of 5 degrees around the electron direction (but ex- 
cluding the radiative one) are included. To suppress the 
background D* + K?e*v., the radiative photon is fur- 
ther required to have a lateral moment within the 
range (0.0, 0.3). This lateral moment, which describes 
the shape of electromagnetic showers, is found in MC 
event studies to peak around 0.15 for photons but to 
vary broadly from 0 to 0.85 for K? candidates. 

In the selection of DT events, the undetected neutrino 
is inferred by studying the missing energy, Emiss, and 
missing momentum, Piss, Which are defined as 


E miss = Ebeam E E, E Ee, (4) 


and 


Pmiss = — [py + pe + pst \/ I seal O° n mp e (5) 


in the rest frame of e*e^ system. Here, E, (Ee) and 
Py (pe) are the energy and momentum of the radiative 
photon (electron), respectively, and mp- is the nominal 
mass of the D^ meson [29]. In calculating piss, only the 
direction vector of the ST D^ candidate, pgp, is used; 
the corresponding magnitude of momentum is fixed. The 
variable Umiss is then defined as 

U miss = miss — IB nasa] (6) 
The distribution of Umiss for the surviving DT candidates 
is illustrated in Fig. [2] the D+ — yetve signals should 
peak around zero, as shown with the dotted curve. 

By studying the MC simulation samples, the back- 
ground from the semi-leptonic decay D* — «?etw, is 
found to have a non-trivial shape in Umiss. Therefore, we 
study the D+ — «?e*v, backgrounds exclusively by se- 
lecting a control sample in data with exactly the same se- 
lection criteria for the ST events and electron candidates 
used in the selection of signal events. The 7° candidates 
are reconstructed from two photons with a 1C kinematic 
fit constraining their mass to the 7° nominal value and 
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FIG. 2. (color online) The Umiss distribution. Dots with error 
bars are data, the red solid-line histogram shows the overall 
fit curve, the blue dash-line histogram shows the background 
D+ — m? etu, and the green shaded histogram includes all 
other background. The black dotted line shows the signal 
MC simulation normalized to the branching fraction B(D* + 


ye* ve) = 100 x 10 5. 


having a fit x? < 20. We extract the yield of the control 
sample DT — ret, Ns. by fitting the corresponding 
Umiss distribution. The expected number of background 
D+ — m?e*v, in the selection of signal D+ > yet, 
N°5”, is calculated with 


NE Nim i 
exp _ DT ST 22,7 
NA = ear (7) 


where Ebr ao is the DT efficiency of Dt > metre, 
DT n° 
Dt + yetve for the tag mode i, individually. The val- 
ues of the corresponding efficiencies are summarized in 
Table[IT] We find NZ = 3016 +68 and N^," = 612 14, 
respectively, where the errors are statistical only. 


E is the rate of mis-identifying D* > net, as 


TABLE II. Summaries of the DT efficiencies of D* — e^ ve 
(ebr) and Dt + m'e*w. (el o), and the rates of mis- 
identifying D* > z? etu, as Dt > yetve [eria 
the branching fraction of Kg — z! and z? > yy are not 
included. The uncertainties are MC statistical only. 


), where 


Tag Mode ene (96) rp di (96) a (%) 

Kiron 27.00 £0.11 27.93 £014 5.32 + 0.07 
Kim m «x? 14.28 + 0.08 13.79 0.11 3.05 + 0.05 
Kon7 28.97+0.10 30.23 +0.14 5.87 + 0.07 
Kn 2? 15.62 + 0.08 15.17 +0.11 3.29 + 0.06 
K9n*z m 17.86 +0.09 17.55 +0.12 3.72 + 0.06 
Kt K-41 21.12 + 0.10 22.28+013 4.19 + 0.06 


An extended unbinned maximum likelihood fit is per- 
formed on the final Umiss distribution as shown in Fig. 
The signal shape is derived from the simulated Dt > 
yetv,. events convoluted with a Gaussian function to 
compensate for resolution differences between data and 


MC simulation. The parameters of this Gaussian smear- 
ing function are extracted according to the discrepan- 
cy in resolution between data and MC simulation in the 
control sample D+ — metve, and are fixed in the fit. 
The shape of the background D* > TÜety, is extracted 
from the simulated D+ — aPe*v, sample, and is nor- 
malized to No. For the other background components, 
the shape from the inclusive MC sample (excluding the 
contribution from D* — aq?e*v,) is adopted and the 
yield is determined in the fit. We obtain a signal yield 
of Npr = —21 + 23, and the resulting branching frac- 
tion is B(D* > yetve) = (—2.5 + 2.7) x 107, where 
the uncertainties are statistical only. Since no obvious 
signal is observed, an upper limit at the 9096 C.L. on 
the branching fraction of Dt — «e*v, will be set be- 
low after taking into account the effects of statistical and 
systematic uncertainties. 


IV. SYSTEMATIC UNCERTAINTIES 


The systematic uncertainties in the selection of the ST 
candidates are assumed to largely cancel, with any resid- 
ual effects being negligible. Other systematic uncertain- 
ties, related to the detection efficiencies, are summarized 
in Table To evaluate the systematic uncertainty re- 
lated to the model of the decay dynamics, an alternative 
signal MC sample based on the single pole model [Bl 
is produced, and the resultant difference in the detection 
efficiency with respect to the nominal value, 3.5%, is as- 
signed as the systematic uncertainty. The uncertainties 
of electron tracking and PID are estimated to be 0.5% 
and 0.596, respectively, by studying a control sample of 
radiative Bhabha scattering events. The uncertainty in 
photon reconstruction is assigned as 1.096, based on a 
study of double-tagged D? — Kgn” events [27]. The un- 
certainty related with the lateral moment requirement for 
the photon is estimated to be 4.4% by studying a photon 
control sample from radiative Bhabha scattering events. 
The quadratic sum of the above systematic uncertainties, 
related to detection efficiency, is 5.896. 

The systematic uncertainty associated with the esti- 
mated number of background D+ — «Pe*r, events in- 
cludes a statistical uncertainty on the size of the DT con- 
trol sample (D* — n? e*v,) of 2.396, and relative uncer- 
tainties on the detection efficiency relative to signal, of 
1.0% for the 7° 1C kinematic fit, and 1.0% for the extra 
photon with respect to the signal. Adding in quadrature, 
the total uncertainty of the background Dt > 7°et1, 
rate is 2.796. Note this value is not the direct fractional 
change in the branching fraction of Dt — yetve, it is 
the fluctuation of background Dt > «?e*v, and will be 
considered along with other effects from the fit procedure. 

Various sources of systematic uncertainties in the 
fit procedure are considered: (a) fits are redone 
with the fitting range being as (—0.15,0.25) GeV or 
(—0.20, 0.25) GeV; (b) the mean and width of the smear- 
ing Gaussian function for the signal shape are varied ac- 
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cording to the corresponding uncertainties obtained from 
the control sample D+ — 7°etv,; (c) the number of the 
background D+ — m?e*vw, is varied according its un- 
certainty (2.796); (d) the shape derived from the inclu- 
sive MC sample is replaced by a second order polynomi- 
al function to describe the other backgrounds excluding 
D+ + nmetve. All of these fitting procedure effects are 
accounted for within the upper limit evaluation described 
next. 


TABLE III. Systematic uncertainties related to detection ef- 
ficiencies in the branching fraction measurement. 


Source Relative uncertainty (96) 
Signal MC model 3.5 
e™ tracking 0.5 
er PID 0.5 
y reconstruction 1.0 
Lateral moment 4.4 
T?) etu backgrounds 2: 


a Note, this value is a fractional change in the m°e+ve rate, not 
in the branching fraction of D+ — yetve. 


V. THE UPPER LIMIT ON BRANCHING 
FRACTION 


'To set the upper limit on the decay branching fraction 
B(D* — yetu), we follow the method in Refs. [271 
which takes into account the effects of both systematic 
and statistical uncertainties. We obtain a smooth proba- 
bility density function (PDF) from the data sample using 
the kernel estimation method [29]. A large number of toy 
MC samples are generated according to the smooth PDF, 
while the number of events in each MC sample is allowed 
to fluctuate with a Poisson distribution according to the 
yield found in the fit to the data sample. The same fit 
procedure used for data is applied to each toy MC sam- 
ple, while randomly making systematic variations in the 
fit procedure, as described in the previous section. In the 
calculation of the branching fraction B(D* — ye ue) for 
the toy MC sample, the DT efficiencies are varied ran- 
domly according to the detection efficiency uncertainties 
(5.896), and the ST yields and the corresponding effi- 
ciencies are varied randomly according to the statistical 
uncertainty due to the size of data and MC samples. The 
resultant distribution of B(D* — ye ue) for all toy MC 
samples is shown in FigB] By integrating up to 90% of 
the area in the physical region B(D* > yet v.) > 0, we 
obtain an upper limit at the 9096 C.L. for the branching 
fraction as B(D* > yetve) < 3.0 x 1075. 


VI. SUMMARY 


In summary, we present the first search for the radia- 
tive leptonic decay Dt — yetve in the charm sector 
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FIG. 3. Distribution of the accumulated branching fraction 
based on toy MC samples generated according to the data. 
The shaded region represents 90% of the physical region. 


based on a DT method using a data sample of 2.93 fb! 
collected with the BESIII detector at a center-of-mass en- 
ergy Vs = 3.773 GeV. No significant DT — yeT us signal 
is observed. With a 10 MeV cutoff on the radiative pho- 
ton energy, the upper limit of the decay branching frac- 
tion for D+ > yetu, is B(D* — yetve) < 3.0 x 10? 
at the 9096 C.L. The result approaches the theoretical 
predictions in Refs. [12]; more data may help to dis- 
criminate among the full suite of theoretical models. 
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